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Abstract 

In  this  paper,  design  and  exergy  method  optimization  of  a  5  kW  power  output  polymer  electrolyte  fuel  cell  (PEFC)  fuel  cell  with  cogeneration 
application  has  been  investigated.  It  is  assumed  that  cooling  water  is  passed  through  the  fuel  cell  cooling  channel,  warmed  up  and  supplied  for 
space  heating  applications.  The  performance  of  the  proposed  system  has  been  optimized  by  the  second  law  based  on  exergy  analysis.  Results 
show  that  for  maximum  system  efficiency  and  minimum  entropy  generation,  fuel  cell  temperature  and  voltage  should  be  as  high  as  possible  in 
the  range  of  application.  Also,  the  fuel  cell  pressure  and  stoichiometric  air:fuel  ratio  should  be  as  low  as  possible  in  the  range  of  application. 
©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  are  electrochemical  devices  that  directly  con¬ 
vert  the  chemical  energy  of  a  reaction  into  electrical  energy. 
The  basic  structure  or  building  block  of  a  fuel  cell  consists 
of  an  electrolyte  layer  in  contact  with  a  porous  anode  and 
cathode  on  each  side.  One  of  the  most  popular  types  of 
fuel  cells  is  the  polymer  electrolyte  fuel  cell  (PEFC).  The 
electrolyte  in  this  fuel  cell  is  an  ion  exchange  membrane  in 
which  corrosion  problems  are  minimal.  A  schematic  repre¬ 
sentation  of  a  PEFC  with  the  reactant/product  gases  and  the 
ion  conduction  flow  directions  through  the  cell  is  shown  in 
Fig.  1. 

From  point  of  view  of  thermodynamic  design,  there  are 
several  methods  to  increase  the  efficiency  of  a  fuel  cell.  The 
most  popular  way  using  power  and  heat  is  combined  cycles 
(CHP)  [1,2].  The  polymer  electrolyte  fuel  cell  which  has  a 
low  power  and  heat  generation  capacity  may  be  used  in  build¬ 
ings  for  local  power  distribution.  For  CHP  application  of  this 
fuel  cell,  water,  which  is  a  cooling  media  in  the  fuel  cell,  is 
used  for  space  heating  or  water  heating  [3].  In  these  fuel  cells, 
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water  is  passed  through  the  fuel  cell  cooling  channel,  warmed 
up  and  is  used  for  the  CHP  application.  There  are  previous 
publications  about  energy  and  exergy  analyses  of  PEFC  fuel 
cell  systems.  These  analyses  have  been  published  with  vary¬ 
ing  degrees  of  cogeneration  in  order  to  seek  and  develop  fuel 
cells  with  better  economic  and  energy  saving  characteris¬ 
tics  than  conventional  power  generating  plants  [4,5].  Several 
hybrid  system  configurations  have  been  suggested  for  effi¬ 
ciency  improvement.  In  addition,  performance  predictions 
under  part-load  conditions  have  been  made  for  the  systems, 
for  which  power  output  is  over  200  kW  [6-8].  Magistri  et  al. 
[9]  proposed  a  very  small  5  kW  gas  turbine  plant  combined 
with  a  SOFC  3 1  kW,  and  evaluated  its  part-load  performance. 
However,  none  of  the  previous  works  have  done  a  complete 
exergy  analysis  of  a  CHP  system,  which  considers  the  heat 
transfer  and  pressure  drop  in  the  reactant,  product  and  cooling 
channel. 


2.  Mathematical  model 

The  proposed  system  under  consideration  is  shown  in 
Fig.  1 .  The  model  consists  of  a  PEFC  fuel  cell  having  5  kW 
power  output  with  cogeneration  application.  In  this  system. 
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Nomenclature 

a 

air:fuel  stoichiometric  ratio 

Cp 

average  specific  heat  (kJ  (kg  K)- 1 

Dh 

hydraulic  diameter  (m) 

e 

exergy  (kJkg-1) 

f 

friction  coefficient 

F 

Faraday  constant 

h 

convection  heat  transfer  coefficient 

(W  (m2  °C)-1) 

iif 

enthalpy  of  formation  (kJ  kg-1) 

k 

heat  conduction  coefficient  (W  (m°C)-1) 

K 

correction  factor  (Eq.  (16)) 

L 

length  of  passage  (m) 

m 

mass  flow  rate  (kg  s- 1 ) 

M 

atomic  number 

n 

number  of  cells 

Nu 

Nusselt  number 

P 

pressure  (kPa) 

A  P 

pressure  loss  (kPa) 

Pc 

perimeter  of  cooling  channel  (m) 

Q 

heat  production  (kW) 

R 

universal  gas  constant  (kJ  (krnol  K)-1) 

Re 

Reynolds  number 

T 

temperature  (K) 

V 

velocity  of  fluid  in  channel  (ms-1) 

V 

fuel  cell  voltage  (V) 

w 

power  (kW) 

X 

mole  fraction 

Subscripts 

ch 

chemical 

f 

fuel  cell 

fi 

inlet  water  temperature 

fo 

outlet  water  temperature 

in 

inlet 

K 

local 

1 

friction 

0 

standard  condition 

out 

outlet 

ph 

physical 

s 

diffusion  layer 

tot 

total 

w 

water 

Greek  symbols 

P 

viscosity  coefficient  (Pa  s) 

P 

density  (kgm-3) 

V 

correction  factor 

P 

system  efficiency 

Anode  Outlet  PEM  Cathode  inlet 


Fig.  1 .  Schematic  of  a  PEFC  fuel  cell  with  proton  exchange  membranes. 


cooling  water  is  passed  through  a  cooling  channel,  warmed 
up  and  supplied  for  space  heating.  The  following  reactions 
take  place  in  the  anode  and  cathode  of  the  proposed  fuel  cell: 

H2->  2H++2e~  (1) 

2H+  +  2e~  +  i02  H20  (2) 

Based  on  the  control  volume,  which  is  shown  in  Fig.  2,  and 
the  mass  flow  rates  of  inlet  hydrogen,  the  inlet  and  outlet  air 
and  steam  can  be  calculated  as  follows: 


mH2,in  -  tFVmh2 

(3) 

W  1 

^  air,  in  —  A  ^.Mair 

4FV  xo2 

(4) 

W  1 

^air,out  —  .  ^ T  l)^air 

4FV  xo2 

(5) 

w 

Wl  steam,  out  —  0 

2r  V 

(6) 

Cooling  water 


Fig.  2.  Proposed  control  volume. 
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Fig.  3.  Flowchart  to  calculate  outlet  temperature  of  the  cooling  water. 


Also,  the  heat  generation  can  be  calculated  by  the  following 
equations: 

Q  -W  =  "»air,out(%  +  CP(7f  -  T0))air  out 

+  ^  '  tn steam, out (hf  +  Cp(7f  —  7o))steamout 

^  '  hrair,in(^f)air,in  ^  '  777 M2, in  ((b) 1 1; .  in  (7) 

It  should  be  noted  that  Eq.  (7)  is  based  on  the  two  assump¬ 
tions,  which  are  given  as  follows: 

1.  inlet  air  and  hydrogen  are  at  standard  condition 
(T0  =  298  K,  P0  =  1  atm); 

2.  outlet  air  and  steam  are  at  fuel  cell  temperature. 


Mass  flow  rate  of  cooling  water  is  calculated  as  follows: 


7«w 


Q 

CP(7f0  -  7fj) 


(8) 


The  unknown  outlet  temperature  of  the  cooling  water  can  be 
calculated  based  on  the  flowchart,  which  is  shown  in  Fig.  3. 
The  specification  of  the  system  is  summarized  in  Table  1.  The 
configuration  of  cooling  channel  is  shown  in  Fig.  4.  Also, 
the  Reynolds  and  Nusselt  numbers  are  calculated  as  follows 


Table  1 

Specification  of  the  system 

W(w) 

5000 

V  (V) 

220 

A  (cm2) 

250 

k 

2 

n 

440 

P(  kPa) 

300 

45(cm) 


50(cm) 


0.5 


mm) 


0.5(mm) 


Fig.  4.  The  configuration  of  the  cooling,  reactant  and  product  channels. 


[10,11]: 

Re  =  4mw  ,  Re  <  2300  (9) 

jrDhun 

Nu  =  7.53  (10) 


The  convection  heat  transfer  coefficient  h  and  outlet  tempera¬ 
ture  of  the  cooling  water  7fo  are  calculated  [12],  respectively: 

k 

h  =  — Nu  (11) 

Dh 


Tfo  =  T{-  (7f  -  rfi)  exp  (  .  PfL  //)  (12) 

Frictional  pressure  loss  in  the  cooling,  reactant  and  product 
channel  of  the  fuel  cell  is  calculated  as  follows  [10,1 1]: 

L  v 2 

A  Pl  =  f—p-  (13) 

Dh  l 

Also,  friction  coefficient  in  the  cooling  channel  is  [10,1 1]: 


(14) 


where  this  coefficient  in  reactant  and  product  channel  is 

[10,11]: 


64 

/  =  P7T  (15) 

Re 

The  value  of  correction  factor  ip  is  equal  to  0.85  [10,11]. 

Focal  pressure  loss  in  the  cooling,  reactant  and  product 
channel  of  the  fuel  cell  is  as  follows  [10,1 1]: 


APk  =  KpV-  (16) 

The  value  of  the  multiplier  K  depends  on  the  channel  geom¬ 
etry  and  the  range  of  its  variation  is  1</>T<2  [10,11]. 
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The  total  pressure  loss  in  the  cooling  channel  can  be  cal¬ 
culated  by  the  following  equation  [10,1 1]: 

APt  =  APi  +  APk  +  APS  (17) 


Also,  pressure  loss  in  the  diffusion  layer  can  be  calculated  as 
follows  [13]: 


384000 

ATj  =  - m 

n 


(18) 


Finally,  the  total  pressure  loss  in  the  fuel  cell  is  as  follows 

[13]: 


L  v 1 
A/)  =  f — p  — 
Dh  2 


v 2  384000 

Kp—  H - m 

2  n 


(19) 


2.1.  Exergy  analysis 


Physical  exergy  is  associated  with  the  temperature  and 
pressure  of  the  reactants  and  the  products  in  the  fuel  cell 
system.  The  physical  exergy  is  expressed  in  terms  of  the  of 
enthalpy  and  entropy  difference  from  those  standard  condi¬ 
tions  of  temperature  and  pressure  T0  =  298  K  and  PQ  =  1  atm, 
respectively.  The  physical  exergy  of  an  ideal  gas  with  con¬ 
stant  specific  heat  C p  can  be  written  as  follows: 


eph  =  CpT0 


+  RT0  In 


(20) 


Tf(K) 


Fig.  5.  Variation  of  entropy  generation  with  fuel  cell  temperature  at  different 
stoichiometric  fuel  air  ratios. 


Assuming  a  constant  number  of  cells  be  equal  to  440  and 
output  power  of  5  kW.  The  entropy  generation  function  can 
be  more  simplified  such  as: 

Sgen  =  f(V,  a,  Tf,  P)  (27) 

3.  Results  and  discussions 


The  chemical  exergy  is  associated  with  the  departure  of  the 
chemical  composition  of  a  system  from  that  of  the  environ¬ 
ment  can  be  calculated  by  equation  below: 

eCh  =  V  v«cch  +  Pin  V  xn  In  xn  (21) 

The  total  exergy  flow  can  be  calculated: 

etot  —  £ph  T  ech  (22) 

As  a  result,  the  entropy  generation  of  the  system  can  be  cal¬ 
culated  as  follows: 


Effect  of  the  fuel  cell  temperature  and  stoichiometric  ratio 
of  fuel  air  on  the  entropy  generation  and  system  efficiency 
are  shown  in  Figs.  5  and  6. 

It  is  obvious  from  the  above-mentioned  figures  that  the 
increase  in  fuel  cell  temperature  causes  the  increase  in  outlet 
temperature  of  the  cooling  water,  which  is  supplied  for  co¬ 
generation  in  this  system  (Eqs.  (11)  and  (12)),  so  that  based 
on  this  effect  the  system  efficiency  increases  and  entropy 
generation  decreases. 

The  increase  in  stoichiometric  air:fuel  ratio  causes  the 
increase  in  outlet  air  mass  flow  rate,  which  is  warmed  up 


Sgen —  y.  |  'y  ’hmCtot.in  ^  '  Woutgtot,out  W  J  (23) 

The  second  law  efficiency  of  the  system  is  as  follows: 

_  mwCp,w(7fo  -  Tfj)  +  W  /0/n 

V  —  ^  :  (24) 

/  M.’^inCini  in  /  .n,,tmout^tot.out 

Considering  Eqs.  (3)-(23),  the  functional  relationship  of  the 
entropy  generation  of  the  system  can  be  written  as  follows: 

•Sgen  =  f(W,  V,  a,  Tf,  T],  L,  pc,  P)  (25) 


Noting  that  the  cooling  and  reactant  and  product  channels 
have  a  constant  geometry  which  is  shown  in  Fig.  4,  the  func¬ 
tion  of  the  entropy  generation  can  be  simplified  as  below: 


Sgen  =  f(W,  V,  a ,  Tf,  P,  17) 


(26) 


Fig.  6.  Variation  of  system  efficiency  with  fuel  cell  temperature  at  different 
stoichiometric  fuel  air  ratios. 
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Tf(K) 


Fig.  7.  Variation  of  entropy  generation  with  fuel  cell  temperature  at  different 
voltage. 

in  the  fuel  cell  and  is  rejected  to  environment.  Due  to  this 
rejection  of  heat,  entropy  generation  increases  and  system 
efficiency  decreases. 

As  seen  in  Figs.  5  and  6,  fuel  cell  temperature  has  a  greater 
effect  on  the  system  efficiency  than  entropy  generation.  Now, 
consider  ainfuel  ratio  a  =  1 ;  the  increase  in  fuel  cell  tempera¬ 
ture  from  330  K  to  550  K  decreases  entropy  generation  about 
21%  but  increases  system  efficiency  about  79%.  Variation 
of  entropy  generation  and  efficiency  of  the  system  with  fuel 
cell  voltage  and  temperature  are  shown  in  Figs.  7  and  8,  re¬ 
spectively.  The  increase  in  fuel  cell  voltage  from  140  V  to 
220  V  at  constant  temperature  decreases  entropy  generation 
about  71%  and  increases  efficiency  about  44%.  This  variation 
shows  the  vital  influence  of  fuel  cell  voltage  on  the  system 
entropy  generation  and  efficiency.  Considering  Eqs.  (3)-(6), 
it  is  obvious  that  the  increase  in  voltage,  decreases  mass  flow 
rate.  This  reduction  decreases  pressure  loss  according  to  Eqs. 


Fig.  8.  Variation  of  system  efficiency  generation  with  fuel  cell  temperature 
at  different  voltage. 


Fig.  9.  Variation  of  entropy  generation  with  fuel  cell  temperature  at  different 
pressure. 

(13)— (19),  whereas  this  variation  decreases  entropy  genera¬ 
tion  and  increases  system  efficiency. 

Figs.  9  and  10  show  the  effect  of  fuel  cell  pressure  and 
temperature  on  the  entropy  generation  of  the  system  and 
efficiency,  respectively.  As  shown,  the  increase  in  fuel  cell 
pressure  increases  entropy  generation  and  decreases  system 
efficiency.  Focussing  on  Eq.  (20),  it  is  understandable  that  the 
increase  in  pressure  increases  physical  exergy  at  the  inlet  and 
outlet  of  the  fuel  cell,  but  this  effect  is  greater  at  the  fuel  cell 
inlet  than  the  outlet,  which  causes  increasing  net  physical  ex¬ 
ergy  in  the  fuel  cell.  According  to  Eqs.  (23)  and  (24)  entropy 
generation  increases  and  system  efficiency  decreases. 

The  most  effective  parameter  affecting  entropy  generation 
and  system  efficiency  according  to  Figs.  5-10  is  voltage  and 
the  pressure  has  a  lesser  effect  on  system  performance,  so  for 
the  maximum  efficiency  and  minimum  entropy  generation, 
the  fuel  cell  temperature  and  voltage  should  be  as  high  as 


Fig.  10.  Variation  of  system  efficiency  with  fuel  cell  temperature  at  different 
pressure. 
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possible;  meanwhile,  the  pressure  and  stoichiometric  airifuel 
ratio  should  be  as  low  as  possible. 

4.  Conclusions 

A  PEFC  fuel  cell  with  a  5  kW  power  output  and  a  cogener¬ 
ation  application  has  been  thermodynamically  designed  and 
its  optimized  performance  has  been  considered  by  exergy 
analysis.  The  effect  of  variation  in  fuel  cell  temperature,  sto¬ 
ichiometric  air:fuel  ratio,  voltage  and  pressure  on  the  sys¬ 
tem  has  been  investigated.  The  following  conclusions  are  ob¬ 
tained: 

1 .  Exergy  analysis  is  a  suitable  tool  for  system  optimization. 

2.  Increase  in  fuel  cell  temperature  and  voltage,  decrease  en¬ 
tropy  generation,  and  as  a  result  the  fuel  cell  efficiency  in¬ 
creases,  whereas  an  increase  in  the  stoichiometric  air:fuel 
ratio  and  fuel  cell  pressure,  increase  entropy  generation, 
and  as  a  result,  the  system  efficiency  decreases. 

3.  The  fuel  cell  voltage  has  the  most  influence  on  entropy 
generation  and  system  efficiency,  whereas,  fuel  cell  pres¬ 
sure  has  the  weakest  effect  on  entropy  generation  and 
system  efficiency.  So,  fuel  cell  temperature  and  voltage 
should  be  as  high  as  possible.  Also,  the  fuel  cell  pres¬ 
sure  and  stoichiometric  air:fuel  ratio  should  be  as  low  as 
possible. 
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